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SUMMARY 


The  need  for  better  dead-reckoning  data  from  several  trials  on 
HMAS  KIMBLA,  which  is  not  fitted  with  a speed  log,  led  to  this  attempt 
tc  model  the  effect  of  wind  and  sea-state  on  ship’s  speed.  Weather 
data  comprised  only  normal  ship’s  log  observations,  which  Imposed  a 
constraint  on  the  factors  considered.  In  particular,  no  specific  data 
on  wave  and  swell  periods  were  recorded.  Equations  were  set  up  to  allow 
for  augmentation  of  resistance  due  to  waves,  swell  and  applied  rudder, 
in  addition,  the  effects  of  the  direct  push  of  wind  and  waves  on  the 
forward  and  trarsvc-rse  components  of  ship's  velocity  were  considered. 

The  parameters  were  fitted  to  observed  data  on  the  set  of  the  ship  in 
varying  weather  conditions.  Some  recommendations  for  future  work  are 
included. 
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1 . INTRODUCTION 

The  effect  cf  wind,  swell  and  waves  on  the  speed  of  a ship  thiough 
the  water  is  quite  complex  and  the  simple  study  presented  here  was  made 
primarily  to  reduce  the  navigational  data  from  two  trials.  The  data  for 
this  particular  study  were  obtained  from  tMAS  KIMBTA  Trial  19/74  during  the  period 
8-13  December  1974.  HMAS  KIMBLA  is  not  fitted  with  either  a distance  or 

speed  log.  If  she  had  been,  it  would  not  have  been  necessary  to  make  the 
study  and  until  she  is  so  fitted,  insufficient  actual  data  are  available  to 
warrant  further  refinement  in  the  theoretical  model.  At  the  time,  however, 
the  need  for  some  current  measurements  in  the  neighbourhood  of  an  eddy 
warranted  a reasonable  attempt  to  deduce  ships'  speed  from  engine  revolutions 
and  standard  weather  observations.  The  principal  reference  for  an  appreciation 
of  the  problem  has  been  the  work  by  Kent  (1958).  His  work  discusses  the 
problem  with  reference  to  ships  somewhat  larger  than  HMAS  KIMBLA,  but 
nevertheless  is  most  instructive.  His  Figure  75  is  shown  as  our  Figure  1, 
illustrative  of  what  may  happen  to  the  speed  of  a low  powered  ship  when  the 
weather  shifts  its  direction  to  her  course  and  about  which  he  makes  the 
following  comments: 

(1)  "In  all  winds  below  50  knots  the  greatest  loss  in  speed  occurred 
with  the  weather  about  30°  on  the  bow,  although  the  differences  in  ship  speed 
were  small  between  'head  on'  and  40°  on  the  bow. 

(2)  For  wind  velocities  cf  50  knots  and  upwards,  the  ship  lost  speed 
rapidly  as  the  weather  direction  approached  'broadside  on'.  This  was 
primarily  due  to  the  rapid  increase  in  helm  angle  necessary  to  keep  the  vessel 
on  her  course,  which  makes  the  rudder  resistance  a major  retarding  force. 

(3)  In  following  wind  and  waves,  the  increase  in  ship  speed  over  that 

in  smooth  water  was  greatest  at  the  lower  wind  speed  of  30  knots,  with  the 
weather  between  'tail  on'  and  30  on  the  quarter.  This  was  due  to  the 

greater  'push'  on  the  stem  by  following  waves  when  they  are  short  by  comparison 
with  the  ship's  length,  than  the  very  much  diminished  push  of  longer  waves. 
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(4)  The  ship  speed  in  30-knot  winds  would  have  been  from  0.3  to  0.5 
of  a knot  more  in  wave  directions  ahead  or  on  the  bow  but  for  the  drop  in 
propulsive  efficiency  due  to  the  greater  augment  in  resistance  in  short  waves  . 
In  higher  winds  the  ship's  speed  would  have  been  less  by  app.oxiraat^ly  the 
same  amount  but  for  the  drop  in  the  augment  of  rasistance  in  w£ve.s  longei.'  than 
the  vesse]  , which  reduced  the  thrust  needed  to  d’^ive  tha  ship  and  raised  her 
propulsive  efficiency. 

(5)  With  the  weather  ahead  or  on  the  bow  the  average  percentages  of 
the  total  resistance  due  to  wind,  hull  and  rudder  wore  as  follows: 


Percentages  of  the  total  ship  resistance  with  the  weather  'head  cn ' . 


Windspeed  (knots) 

30 

40 

50 

Hull  resistance  in  rmooth  water 

28 

21 

Eytra  hull  resistance  in  rough  water 

36 

34 

Wind  resistance 

10 

16 

24 

32 

39 

Rudder  resistance 

1 

2 

L_ 

4 

6 

As  u.a  weather  direction  veered  towards  the  beam,  the  percentages 
of  the  wind  and  the  extra  hull  resistances  dropped  sharply  and  that  due  to 
the  rudder  rose  rapidly". 

He  then  discusses  the  loss  of  ship  speed  with  change  in  wave  length, 
something  which  our  model  could  not  pursue  owing  to  lack  of  recorded  data. 

The  approach  used  here  was  to  write  down  a set  of  equations 
involving  only  the  available  weather  data  in  order  to  compute  the  modifying 
effact  of  the  weather  on  the  smooth  water  speed  of  the  ship.  It  was  then 
hoped  that  a critical  appraisal  of  all  the  data  (bridge  log,  fixes,  expected 
currents  etc  from  KIMBLA  Trial  19/74)  would  enable  some  actual  values  to  be 
assigned  to  the  constants  in  the  equations.  This  was  achieved  by  reasonable 
estimates  and  trial  and  error,  but  a more  sophisticated  multivariate  correlation 
analysis  could  be  applied  if  the  need  and  the  data  (as  well  as  the  equations) 
warranted  it. 

The  equations  were  set  down  with  some  cognizance  of  the  physical 
forces  Involved  - e.g.  one  might  expect  that  the  energy  transmitted  to  the  ship 
head  cn  into  raves  would  increase  as  the  square  of  the  wave  height  - but  they 
are  not  meant  in  any  way  to  be  regarded  as  analytic  solutions  to  the  overall 
problem.  The  approach  is  largely  empirical.  it  became  apparent  to  the 
author  when  writing  up  this  work  that  the  analysis  and  model  calculations  could 
be  improved  in  several  ways.  These  have  been  noted  appropriately  in  the  text. 
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2.  FORM  OF  EQUATIONS 

The  ustial  pertinent  weather  observations  consist  of  wind  speed  .ind 
direction,  sea  height  and  swell  height  and  dir^iction.  Sea  and  swell  wave- 
length is  not  noted,  so  that  possible  resonance  effects  cannot  be  considered. 

Use  the  following  notation  and  units: 

^ = actual  ■ '.p's  velocity  through  the  water  (kn) 

* ship's  speed  in  absolutely  calm  weather  (kn; 
h ® observed  wt  -'e  height  (m) 

V. 

h “ observed  swell  height  (m) 

8 

w wi  id  speed  (kn) 

= diI.':.t^  1 from  which  wind  and  "sea  are  coming  (degrees  true) 

, from  which  .swell  is  coming  (degrees  true) 

* ct  . Steered  by  the  ship  (degrees  true) 

= cci.fce  made  good  through  the  water  (degrees  true) 

In  addifioQ,  resolve  ^ into  two  components,  vj  in  the  direction  i|)*  and  a 
normal  co..sponent  V2  • 


Write 

0'  => 

w 

1 

> 

o 

and 

6*  » 
s 

0 - 
s 

as  -‘he  relative  directions  (to  the  ship's  bow)  of  wind  and  swell  respectively, 
such  that 

I 0’  i S 180°. 


Assume  that  the  effects  of  the  weather  can  be  treated  separably  in  the  form 


Vl 

* K ®S  ’ ®W» 

O W S ^ w 

....(1) 

V2 

CM 

• 

• 

« 

• 

where  the  function  f can  be  written 

• ^2(w,  e’  ) . f3(h  , 6’  ) ....(3) 

W W W 8 o 


f 


io 


4 


t 

i 

! 


Thus,  we  have  defined  retardation  (or  enhancement)  factors  f^,  f2  and  f3  for 
wave,  wind  and  swell  effects  respectively,  acting  on  the  forward  speed  v^ 
and  rn  additional  additive  function  £4,  due  to  wind  and  dependent 
on  ship’s  velocity,  also  contributing  to  vj  . The  normal  compcnent  of  ship's 
speed  V-  (the  leeway)  is  defined  by  £5  as  a function  of  vj  and  wind  velocity. 

The  question  now  arises  as  to  the  form  of  the  functions  fj  to  £5.  Fa^'tors 

fj,  to  and  £3  form  a similar  group,  so  let  us  consider  the''  together.  A.ssume 
that  the  power  delivered  from  the  ship’s  engine  is  dependent  only  on  the  product 
of  engine  revoluttons  and  shaft  torque  and  is  independent  of  variation*,  in  ship’.s 
speed  (HMA£  KKiBLA,  for  example,  .sets  a constant  steam  pressure  and 

throttle  setting  is  only  changed  when  different  revolutions  are  called  for.) 

The  effective  power  i-s  used  to  drive  the  ship  forward  at  speed  vj  and,  in 
general, 

P = F V, 
e ^ 

whev'o  F is  the  forward  thrust.  Hence,  for  constant  fhcottie  setting 

F ® --  for  weather  no  worse  than  roodorate. 
vi 

As  the  sea  state  increases,  Vj  decreases  and  F increases.  This  corres- 
ponds to  an  increasing  propeller  slip  ratio  and  increasing  angle  of  Incidence 
of  the  water  on  the  propeller  blade;  so  F increases.  (See  Todd,  1967). 

Now  consider  the  equilibrium  between  forward  thrust  F and  the 
overall  drag  or  resistance  D of  the  ship  to  forward  action.  Write  D in 
the  foim 

D « A vi^ 

so  that  C.  takes  the  form  of  a drag  coefficient,  A is  son-c  const '.nt 
d 

(dimension  i-)  for  the  ship. 


Equating  F and  D gives  us 


constant 


-1, 


73 


4 


J 


or 


....(4) 
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Call  Cj  the  resistance  coefficient  which  will  vary  with  the  weather  (wind, 
sea  state  etc.).  Suppose  in  calm  weather  it  has  the  value 


Write  C,  = (1  + 3)  (5) 

d do 

..  •#.  •••  • • •*  •• 

where  0 represents  an  augmentation  of  our  resistance  coefficient  due  to 
weather.  Thus  in  calm  weather  0=0  and  vj  takes  a value  say. 

Substituting  in  (4)  we  obtain 

-Va 

vi  = V (1  + 0)  .... (6) 

o 

This  is  the  form  in  which  we  shall  write  the  retardation  factors  fj,  £2  and  fs. 
3.  RETARDATION  DUE  TO  WAVE  RESISTANCE 


The  function  fi(h  , 6'  ) is  assumed  to  have  the  form 
‘ w w 


f 1 = 


...(7) 


where  is  an  empirically  set  normalization  constant  andg^(6’^)  describes 

the  variation  with  incident  angle,  .such  that  g(0)  = 1.0  . 

(The  value  of  the  index  used  in  '"/)  in  the  model  calculations  was  actually 
-O.IS  instead  of  -0.33,  but  the  difference  is  not  significant.) 

When  h = H and  6'  =0,  fi  = 0.79.  H tras  finally  set  at  2.44  m 

WWW*  w 

(8  ft). 

The  function  g^(6|^)  controlling  the  variation  of  fj  with  relative 
angle  of  incidence  is  shown  in  Figure  2,  along  with  the  corresponding  function 
for  swell  effect,  g (6'),  in  the  factor  13.  These  functions  have  been  set  up 
entirely  on  the  basis  of  "reasonableness",  lacking  much  in  the  way  of  observed 
data.  Note  that  the  effect  of  swell  is  assumed  to  be  zero  once  the  incidence 
is  abaft  the  beam.  Some  retardation  from  wave  action  is  still  assumed,  how- 
ever, for  1 9^1  > 90°.  This  will  be  counteracted  to  a varying  degree  by  the 
positive  push  allowed  for  in  the  function  fi,,  augmentation  of  resistance  and 
push  being  considered  independently.  The  differences  between  swell  and  wave 
action  are  introduced  because  of  the  difference  in  wavelengths. 


i 
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A.  RETARDATION  CAUSED  BY  AFPI.IED  RUDDER 


According  to  Kent,  under  strong  beam  winds,  the  applied  rudder 
necessary  to  hold  a given  course  can  cause  sufficient  resistance  t^  dominate 
the  other  resistance  effects.  It  may  not  even  be  possible  to  held  a course 
'“in  beam  winds  of  above,  say,  50  knots,  ad' illustrated  in  Figure  1'.  Such 
retardation  is  shown  as  reaching  a maximum  somewhat  abaft  the  beam,  at  about 
110°  from  the  bow.  This  would  also  accord  with  the  more  usual  situation  in 
which  steering  is  more  difficult  with  the  seas  coming  from  abaft  rather  than 
forward  of  the  beam.  Greater  use  of  the  rudder  will  increase  the  average 

resistance. 


Hence,  rather  crudely,  we  define 


0'  = • |e‘  I for  1 0’  I < 110^ 

r 11  ' w ' ' w ' 


...(8a) 


0’  = |0'  I - 20°  for  0’  5 110° 


(8b) 


so  that  sin  0'  = 0 for  0’  = 0 , 

r w ’ 

= 1.0  for  I 0'  I = 110°  (maximum  effr  .t). 
' w' 

= 0J4  for  I = 180°  (wind  astern) 


and  now  define  f2(w,  0'  ) by 


f2  = 1 + C^(w  sin  e’j.  )^  ] 


where  is  an  empirically  set  constant. 


....(9) 


The  terra  w sin  0^  was  raised  to  the  third  power  i iside  the  bracket 

as  a step  towards  satisfying  the  general  form  of  the  crfpi  t in  Figure  1. 

Note  the  rapid  Increase  in  this  effect  shown  by  Kent  as  the  wind  inc eases  from 

40  to  50  knots.  Also,  the  effect  is  strong  with  0*  between  80°  and  140°  . 

w 

The  whole  term  was  Included  with  the  index  set  at  -0.5  instead  of  -0.33 
(suggested  by  (6))  to  further  strengthen  the  function.  The  resultant  form  of 
f2  as  defined  by  (9)  and  used  in  the  calculations  is  showii  tu  Fig.  3 for 
winds  varying  from  0 to  51  knots.  In  retrospect,  it  might  have  been  better 
to  leave  the  index  at  -0.33  and  to  lise  w and  sin  0 senarately  to  higher 
powers  to  better  match  Figure  1, 
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A consideration  jf  the  data  along  part  of  the  track  of  Trial  \^llk  (leg  V -W, 
see  Figure  6 ) , in  which  beam  winds  of  25  knots  were  encountered,  showed  that 
the  effect  did  not  appear  marked  under  those  conditions.  Herce  was  set 
at  bclO"^  . 


5 . RETARDATION  DUE  TO  SWELL 


This  is  similar  to  factor  Define 


f3 


1 + 83(6;) 


g (O')  is  shown  in  Figure  2. 
s s 


... .(10) 


There  were  insufficient  data  to  really  distinguish  between  Che  effects  of 

waves  and  swell,  but  as  the  ship's  log  records  these  Independently,  witii 

directions  noted  for  each,  the  factors  f^  and  fg  hc.ve  been  separately  included. 

Figure  4a  shows  a plot  of  swell  against  wave  height  observed  during  Trial  19/74. 

Swell  height  seems  roughly  correlated  with  wave  height  such  that  h = 2 h . 

s w 

Hence  the  normalization  height  in  (10),  H , has  been  set  at  twice  the  value  used 

s 

for  in  (7).  Thus  the  factors  fj  and  tend  to  bo  of  the  same  magnitude 
under  average  conditions.  This  is  also  consistent  with  the  observation  that 
waves  (coming  from  bow  on,  say)  have  a greater  retardation  effect  than  swell  of 
the  same  height  because  the  former  have  shorter  wavelengths  and,  being  steeper, 
lose  more  energy  to  the  ship.  Also,  the  separation  of  the  sea  state  into  swell 
and  waves  is  a matter  of  visual  judgement  by  the  ship’s  officers  and  is  often 
most  imprecise.  Insofar  as  Trial  19/74  is  concerned,  the  absolute  accuracy 
of  the  estimati'^ns  of  sea  and  wind  state  is  not  important,  as  the  various 
constants  in  the  model  equations  were  adjusted  to  suit;  however,  if  this  model 
were  used  with  data  from  other  trials  the  abso.lute  accuracy  of  the  weather 
observations  would  be  important.  Fortunately,  as  the  weather  observations  are 
made  by  a number  of  watch-keeping  officers,  individual  tendencies  to  under  or 
over-estimate  the  parameters  are  probably  evened  out. 


6.  THE  PUSH  OF  WIND  .4ND  WAVES 

Apart  from  the  retardation  already  discussed  which  the  ship  suffers 
due  to  travelling  through  a disturbed  sea,  the  vessel  also  receives  a direct 
push  due.  to  Incident  wind  and  waves.  This  can  be  conveniently  resolved  into 
two  components,  parallel  and  normal  to  the  ship's  heading.  At  the  time  of  this 


study  and  for  the  purpose  of  these  foirmulae,  there  seemed  little 
point  in  trying  to  distinguish  between  wind  and  waves  in  the  open  ocean,  as 
they  are  usually  quite  well  correlated  and  a common  direction  is  recorded  for 
both.  In  addition,  it  is  not  uncommon  in  practice  for  one  to  be  estimated 
from  the  other.  Usually  the  wind  is  estimated  both  in  magritude  and  direction 

from  the  observed  sea.  Thus  the  recorded  wind  should  be  an  absolute  measure- 
ment, independent  of  the  ship's  motion.  Should  the  wind  be  recorded  by 
instrumentation  on  the  ship,  however,  the  appsrent  velocity  would  have  to  be 
corrected  for  ship's  motion.  The  relation  between  recorded  wave  height  and 
wind  speed  for  Trial  19/74  is  shown  in  Figure  4b.  A reasonable  estimate  of 
wave  height  can  be  obtained  from 

h Cmetres)  = 0.06  x w (knots)  ....(11) 

Supposing  that  the  swell  reached  twice  this  figure,  the  heights  are  still 
considerably  less  than  those  expected  of  a fully  developed  sea.  This  is 
consistent  with  the  rapidly  changing  nature  of  the  weather  at  the  time. 

The  net  addition  to  the  ship's  forward  speed  is  assumed  to  be 


f4(w,  0*  ) = p cos  (6’  + 180°) 

w w 


where  p = C w for  w < 25  kn 

P 

or  - (O'  w)°*5  for  w > 25 kn 

P 

and  C C'  are  fitted  constants. 

P P 


(12) 


The  change  in  the  form  of  p(w)  at  25  kn  is  based  on  the  concept  that 
in  strong  winds  the  waves  get  longer  as  well  as  hip''er  and  thus  the  pushing 
effect  of  the  waves  drops  below  a linear  increcso.  A half-power  law  for  all 
wind  speeds,  however,  over-estimates  the  effect  of  light  winds. 

In  retrospect,  the  following  analysis  is  more  soundly  based; 

From  Section  2,  we  have-  for  the  propulsive  force  on  the  ship 


2 


and,  under  given  weather  conditions,  assume  the  drag  coefficient  C,  to  remain 

d 

constant. 

2 

F = C'vj  , where  C'  is  some  constant.  Therefore 
AF  = 2C*vi  A vj 

Suppose  AF,  due  to  wind  and  waves,  is  of  the  form 
AF  = p(h^,  w,  V) 


= q(h  . w').  r(e',) 

w w 


(14a) 


where  w'  is  the  relative  wind  speed.  Function  q can  reasonably 

be  made  dependent  on  w'^,  but  the  dependence  on  waveheight  is  more  uncertain, 
owing  to  thevariatlon  of  wavelength.  Suppose  for  simplicity, 

q(h^,  w')  = 

The  variation  of  is  not  a simple  cosine  as  used  in  (12).  As  the  angle 

of  the  apparent  wind  to  the  bow  0^,  varies  from  zero  the  projected  area  of  the 
ship  Increases  rapidly  and  r(9\)  Increases  to  a maximum  at  about  30°  off 

V 

the  bow.  Such  a function  is  given  by  Todd  (1967)  reporting  the  work  of  Hughes 
(1930).  Using  this  function,  we  have 

P = 9(\.  w’>  . r(0^,) 

and  i4(h„»  w,  0^,  y)  = Avj 


2C’vi 


(14b> 


so  long  as  the  ship  is  under  way.  Nov  suppose  the  engine  is  stopped  in 
moderate  winds  and  AF  contributes  the  total  propulsive  force. 

We  have  AF  = C ’ ( AV^ ) ^ 

hence  Avj  ®(-^)^  (14c) 

0 

(Sven  at  low  speeds  the  Reynolds  number  is  large  and  the  flow  will  be  turbulent.) 
This  shows  why  the  ship  may  move  through  the  water  at  several  knots  with  moderate 
winds  and  seas  astern,  but  no  engine  power;  whereas  the  additive  effect  xmder 
normal  steam-'ng  conditions  is  much  less  and  decrease?  with  Increasing  speed. 

It  rJLso  suggests  a way  of  evaluating  p/C.  With  engines  off,  Av^  could  be 
measured  under  varying  conditions  in  which  (14c)  applies.  Then  (14b)  could 
be  evaluated.  For  example,  suppose  with  25  kn  winds  astern  Av^  » v^  is 
measured  at  3 kn.  The  addition  under  slmllox  conditions  while  under  way  at  9kn 


10 


would  be,  by  {14b),  only  0.5  kn.  A full  and  proper  analysis  of  these  effects 
is  beyond  the  capacity  of  this  paper  (and  of  the  author).  The  calculations 
were  made  using  (13),  however,  before  Eqs.  (14)  were  appreciated.  The  errors 
are  not  so  severe  as  to  warrant  repeating  the  analysis  on  this  account  * 

An  effective  push  while  under  way  of  1 kn  in  head  (or  tail)  winds 

of  25  kn  seemed  reasonable  at  the  time  and  did  not  appear  inconsistent  with 

the  data,  so,  using  (13),  C = C = 0.04  . In  view  of  (14),  this  value  of  fu 

P P 

(while  under  way)  now  appears  a little  large;  however  it  should  be  noted  that 
Kent  (see  Figure  1)  appears  to  predict  larger  values. 

The  normal  (sideways)  component  of  this  effect  has  an  interesting 
aspect.  Elementary  hydrodynamical  considerations  suggest  that,  provided  the 
ship  is  not  stalled  in  the  water,  the  sideways  speed  is  also  reduced  linearly 
as  the  forward  speed  increases.  Thus  under  normal  conditions,  the  leeway 
angle  should  fall  off  as  the  square  of  the  forward  speed.  Consider  Figure  5. 
We  suppose  that  the  ship  is  travelling  through  the  water  with  a velocity 

‘ [v^] 

The  ship  car.  be  considered  stationary  with  the  water  stream  incident  with 
velocity  -V  . The  ship  makes  leeway  angle  a , given  by 

-1  ''2 
a = tan  ^ — 

1 

Suppose  the  sideways  force  due  to  wind  and  sea  is  F^  . Under  equilibrium 
conditions  this  will  bo  .Jdlanced  by  a hydrodynamic  sideways  ’lift*  L,  which  for 
small  a is  given  by 

■ .2 

L « a V ..,.(15) 

I .MHV  » 

F^  depends  on  the  wind  and  sea  state  and  the  relative  heading  of  the  ship. 

For  given  conditions,  F^  (and  hence  L)  is  independent  of  j » so  from  (15) 
we  have 

a « ~^2  •••• 

i vl 

To  a sufficient  approximation,  we  can  substitute 


Now,  for  the  model,  initially  define  the  -ormal  velocity  component  V2  by 


V2  = C p sin(0'  + 180°)  ....(19) 

^ w 

where  p is  defined  as  in  (13)  (or  preferably,  for  new  work,  by  an  equation 
similar  to  (14a)  using  absolute  wind  velocity)  and  is  a constant. 

From  (1)  and  (3)  we  have 

VI  = v^  . fi  . f2  . f3  + U 
I I o 

If  |v2l  > vj  (ie  a > 45  ) the  ship  is  almost  certainly  hydrodynainically 
stalled  and  there  is  no  justification  for  reducing  the  initial  value  V2 

according  to  (18).  Somewhat  arbitrarily,  write 

V2  = V2  for  jv^l  ^ vj  (stalled)  ....(20e) 

or,  < 

V2  ^ 

\’2  = ~ for  |V2i  Vj  (underway)  ....(20b) 

A different  transition  (stall)  angle  could  be  chosen,  or  a non-linear  relation- 
ship between  L and  a (15)  could  be  incorporated,  but  in  the  absence  of 
measured  data  on  leev;ay  angles,  such  refinements  are  not  presently  justified. 

After  some  consideration,  v;e  set  C = 2,5  . This  gives  results  that  do  not 

9 

conflict  with  impressions  gained  over  several  trial..  According  to  this 
value  a 25  kn  wind  broadside  I'u  to  a stopped  ship  would  cause  it  to  drift 
do\^nwind  at  2.5  kn.  This  may  be  a little  high,  but  the  stalled  condition  was 
not  of  prime  importance.  Underway  at  10  kn  however,  with  the  same  wind  on 
the  beam^V2  reduces  to  0.25  kn  and  the  leeway  to  1,4°. 


' • SHIP’S  SPEED  IN  SMOOTH  WATER 

In  the  absence  of  any  data  to  the  contrary,  it  was  assumed  that  the 
apparent  propeller  slip  ratio  was  constant  over  the  norcsal  operating  range  of 
engine  revolutions.  This  is  defined  by 

* ‘-fe  ••••«'> 


where  P is  the  propeller  pitch  and  n the  propeller  revolutions  (in 
appropriate  units).  A constant  value  of  S implies  a simple  linear  relation- 

Qi 

ship  between  v^  and  engine  revolutions,  ie 


K can  be  regarded  as  constant  over  the  course  of  any  one  trial,  but  its 
value  will  vary  with  time  depending  on  the  state  of  the  ship's  hull,  loading, 
state  of  the  propeller  etc.  The  degree  of  hull  fouling  is  usually  the  major 
factor.  HMAS  KIMBLA,  with  a clean  hull,  usually  works  with  a value  K = 15.0 
r.p.m./knot.  Sufficient  data  existed  from  Trial  19/74  to  make  an  independent 
determination  of  this  figure.  This  is  given  in  Section  9. 

' 8.  PROGRESS  THROUGH  THE  WATER 

The  progress  through  the  water  is  given  by 

2 I 

|v|  = (vj^  + V2  ....(23) 

in  the  direction 

(}i  = ())’  + a 

where  a is  the  leeway  angle. 

9.  FITTING  THE  MODEL  TO  THE  DATA 

Figure  6 shows  a plot  of  the  ship's  track  for  Trial  19/74.  A total 
of  66  fixes  were  obtained  from  a satellite  navigation  system  (accuracy  =0.3 
n. miles)  during  the  course  of  the  trial  (114  hr).  The  fixes  of  significance 
to  this  work  are  lettered.  Associated  with  each  labelled  fix  is  a figure 
giving  the  number  of  hours  into  the  trial  starting  from  0001  lirs(time  zone  K) 
on  08  December  1974.  Figure  7 shows  the  dynamic  topography  of  the  area, 
derived  from  Expendable  Bathythermograph  (XBT)  data  taken  during  the  trial. 

The  expected  surface  currents  are  inversely  proportional  to  the  contour 
spacing. 

Now  the  determination  of  the  model  parameters  rested  on  two  propositions 
Firstly,  that  the  surface  currents  ever  the  legs  CDEF  were  small  enough  for  the 
net  effect  to  be  ignored.  Secondly,  that  the  average  surface  current  over  the 
leg  RS  equalled  that  over  the  leg  ST.  The  reasonableness  of  these  propos- 
itions can  be  judged  by  referring  to  Figure  7,  bearing  in  mind  that  a topo- 
graphic gradient  of  10  dyn.  cm.  (0/  900  ra)  per  18.5  km  (10.0  n ‘miles)  is 
estimated  to  produce  a surface  current  of  0.75  m sec~^(1.46  kn)  . This  figure 
includes  a factor  of  1.15  tc  allow  for  the  fact  that  the  dynamic  ocean  structure 


extends  below  900  m,  but  does  not  Include  any  enhancement  due  to  the  curvature 
of  the  topographic  contours. 


Table  1 ; Details  of  weather  during  Trial  19/74 


Time 

w 

w 

R 

w 

0 

s 

R 

s 

hr 

0 

true 

kn 

ft 

0 

true 

ft 

19.0 

200 

8 

1 

180 

4 

23.0 

180 

14 

1 

180 

6 

27.0 

- 

0 

1 

180 

6 

31.0 

040 

5 

2 

180 

35.0 

000 

19 

2 

020 

39.0 

040 

19 

3 

020 

5 

43.0 

000 

30 

4 

030 

7 

47.0 

010 

30 

6 

010 

12 

51.0 

010 

19 

3 

010 

8 

55.0 

180 

14 

4 

020 

6 

95.0 

020 

5 

0 

040 

3 

99.0 

010 

2 

1 

020 

2 

103.0 

240 

8 

3 

- 

0 

107.0 

230 

14 

2 

250 

4 

111.0 

190 

19 

2 

190 

4 

Table  1 shows  the  recorded  weather  during  those  sections  of  the  trial 
pertinent  to  this  work.  Original  units  are  given.  The  weather  at  the  start 
of  the  trial  was  relatively  light  (winds  $ 16  kn)  from  varying  directions 
until  point  E was  reached,  after  which  the  wind  freshened  to  30  kn  from  010° 
around  point  E* , with  sea  and  swell  reaching  1.8  m (6  ft)  and  3.7  m (12  ft) 
respectively  from  the  same  direction.  At  F’  the  wind  quietened  and  came  in 
from  the  south.  During  legs  RST  the  winds  were  again  less  than  16  kn  from 
varying  directions.  Thus,  variations  in  the  parameters  of  the  model  only  had 
a small  effect  on  the  calculated  progress  over  legs  RST,  but  a considerable 
effect  over  the  leg  EF,  particularly  as  the  wind  and  sea  were  almost  head  on 
at  that  time. 

In  view  of  the  above,  an  iterative  procedure  was  adopted  to  fit  the 
various  parameters  to  the  data.  First  estimates  were  made  and  the  set  of  the 

ship  calculated  along  the  various  legs.  If  the  data  and  model  were  perfect, 
the  calculated  set  between  successive  fixes  would  be  due  to  surface  current 
plus  any  net  error  in  the  fixes.  This  calculated  set  has  been  termed  the 
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\! 


'drift'  in  this  discussion.  In  practice,  this  will  also  contain  errors  in 

the  calculated  dead  reckoning  (DR)  data.  The  value  of  K (Eq.(22))was  adjusted 

to  equalize  the  magnitude  of  the  drift  over  the  legs  RS  and  ST.  Row,  with 

respect  to  the  legs  CDEF,  errors  in  calculated  ship's  speed  will  show  up  directly 

in  the  comparison  between  distance  travelled  through  the  water,  D , and  distance 

w 

travelled  over  the  ground,  D , assuming  no  net  effect  due  to  surface  currents. 

8 _ 

The  criterion  for  the  model  was  to  minimise  ) (D  - D ) over  the  three  legs  CD, 

^ w g 

DE  and  EF,  separately  and  totally.  Having  adjusted  the  model  appropriately, 

the  drift  over  RST  was  recalculated,  K readjusted  and  so  on.  This  process 

conve  ged  after  a few  cycles,  leading  to  K - 15.0  1 0.2  r.p.m./kn  - a value 

ident  '.cal  to  that  traditionally  used  by  HMAS  KIMBLA's  officers.  The  drifts 

during  legs  RS  and  ST  are  0.75  m.sec~^  at  287°  and  286°  true  respectively. 

The  dynamic  topography  suggests  a mean  current  of  about  0.8  - 0.9  m.sec"^  at 

285°  true,  so  the  calculated  directional  data  (for  which  there  were  no  direct 

constraints)  are  most  encouraging.  Interestingly,  the  exponents  of  the 

factors  fi  and  fi,  ((7)  and  (10))  were  initially  each  set  at  -0.5,  due  to  a lack 

of  appreciation  at  the  time  of  (6).  The  data  showed  -0.5  was  too  severe  and 

the  magnitude  of  the  exponents  was  reduced  to  -0.35  . This  gave  a better 

fit  and  now  can  be  seen  to  accord  with  the  value  of  - /3  suggested  by  (6). 

The  values  finally  accepted  for  '[(D  - D ) for  legs  CDEF  are  given  in  Table  2, 

w g 

along  with  the  time  averages  of  v^,  f,  f4  and  f5  . 


Table  2 

Functional  values  over  the  section  CF 


(units  are  hr  and  n. miles) 


leg 

C-D 

D-E 

E-F 

IC-F 

79.2 

32.9 

88.8 

200.9 

.969  0.15 

.970  - .12 
0.854  - .84 


25.50 

33.73 

37.10 


Unfortunately,  these  figures  convey  an  impression  which  is  too 
optimistic.  We  should  be  able  to  obtain  similar  agreement  for  any  portion  of 
the  track  along,  say  BDEF*  . There  is  the  limitation,  however,  that  the  weather 
observations  were  only  made  every  4 hr  and  the  weather  changed  quite  rapidly 
over  this  period.  Hence,  one  might  not  expect  agreement  over  any  leg  much 
shorter  than  = 8 hr.  For  this  reason  not  much  notice  was  taken  of  the  short 
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leg  DE  compared  with  the  two  longer  ones  in  adjusting  the  model  and  the  quite 
good  agreement  shown  in  Table  2 for  leg  DE  is  of  little  moment.  As  a test, 
the  section  BDEF’  has  been  subdivided  in  a different  mr.nner  to  that  used  in 
Table  2.  Using  the  same  model,  the  values  for  ^g)  ^nd  the  time  averages 

of  the  main  functions  were  recalculated  for  the  new  divisions  and  are  given  in 
Table  3.  Each  division  occupies  about  8 hr. 


Table  3 


Functional  values  over  the  section  BF ' . 
(Units  are  hr  and  n. miles) 


1 start 

1 time 

leg 

1 D 

)D  - D 1 At 
^ w g ! 

1 <v  > 1 
error  ‘ o | 

< f > 

<f4> 

<f5> 

1 

' 21.80 

B - C 

70.5 

0.8 

7.53 

0.11 

9.89 

.974 

-.17 

0.04 

I 29.33 

C-  E 

75.2 

-1.0 

' 7.77 

-0.13 

9.80 

.965 

0.09 

-.09 

i 37.10 

E - E' 

45.9 

8.9 

. 8.70 

1.02 

8.02  ! 

.884 

-.8] 

0.13 

, 45. 8C 

E'-  F' 

61.6 

-9.0 

8.50 

-1.06 

8.13  i 

.839 

-.65 

0.11 

Ib-  F' 

253.2 

-0.3 

::2.50 

-0.01 

To  satisfy  the  observed  values  of  distance  travelled  D , the  values  of  f 

g 

legs  EE'  and  E'f’  would  have  to  be  0.76  and  0.97  respectively. 


It  is  at  once  apparent  that  although  the  net  value  of  ^(D  - D ) over 

w g 

the  section  EE'F*  is  satisfactorily  small  (-0.1  n. miles),  the  errors  in  each 
leg  separately  are  quite  gross  (8.9  and  -9.0  n, miles  respectively).  These 
correspond  to  errors  in  ship's  speed  of  about  1.0  kn.  To  take  a closer  look 
at  these  errors,  suppose  for  a start  that  the  values  of  f4  (the  component  of 
the  push  received  from  the  wind  and  waves  along  the  ship's  course)  are  correct. 

for 

The  figures 

for  the  observed  weather  however,  indicate  that  the  winds  were  about  the  same 
and  the  seas  (from  nearly  head-on)  were  more  severe  during  the  latter  leg, 
resulting  in  a lower  value  of  f.  (Each  weather  observation  is  assumed  to  be 
a mid-interval  observation,  ie  to  apply  for  2 hr  before  and  after  the  noted 
time.)  Adoption  of  (14)  to  calculate  f4  might  improve  matters  a little,  but 
no  amount  of  sensible  adjustment  of  the  parameters  in  the  model  will  fully 
satisfy  these  data.  A major  error  in  the  fixed  position  of  E'  would  provide 
an  easy  explanation  for  this  discrepancy,  but  an  examination  of  that  and  neigh- 
bouring fixes  gives  no  evidence  of  such  an  error.  With  the  limitation  of  both 

the  model  and  the  data  we  should  rest  with  the  overall  minimization  of  7(0  - D ). 

^ w g 
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The  values  for  the  various  parameters  contained  in  the  model  equations 
were  obtained  from  the  procedure  discussed  in  Section  9.  These  values  have 
been  given  with  the  appropriate  formulae  in  Section  3-7.  Figure  8 shows  the 
calculated  variation  of  speed  v^  for  KMAS  KIMBLA  (in  the  direction  of  ship's 
head)  as  a function  of  wind  velocity  and  relative  heading.  It  is  of  interest 
to  compare  this  with  Figure  1,  which  shows  the  variation  for  a large  vessel 
given  by  Kent  (1958).  The  calculations  for  Fig. 8 were  made  assuming  that 
wind,  waves  and  swell  all  come  frcm  the  same  direction  and,  furthermore,  that 
wave  height  is  proportional  to  wind  speed  (given  by  Eq.(li))  and  swell  height 
is  twice  the  wave  height. 

HMAS  KIMBLA  has  the  following  measurements: 1 122  tonnes  (1002  ions) 
displacement,  45.7  m (i50  ft)  length  overall,-  9.8  m (32.2  ft)  beam,  2.7  in 
(8.7  ft)  draft  forward  and  4.6  m (15,2  ft)  draft  aft.  Power  is  350  I.K.P.  at 
175  r.p.m.,  a speed  not  usually  attained  except  perhaps  within  24  ht  of  reaching 
home  port.  Figure  8 iiidicates  a larger  initial  susceptr.bility  to  sea  sUte 
than  that  shown  in  Figure  1.  A 30  kn  wlnd/sea  head  on  apparently  reduces 
HMAS  KIMBLA's  speed  from  10  kn  (at  150  r.p.m.)  to  about  6.2  Un,  a drop  of  38%j 
whereas  the  reduction  for  the  larger  vessel  is  from  14.4  xi.  to  12.1  kn,  ie  16%. 
This  is  not  surprising,  considering  that  the  larger  ship  has  about  six  times 
as  much  power.  The  most  notable  difference  between  the  two  is  that  of  the 
apparent  behaviour  in  strong  beam  winds.  As  discussed  in  Section  4,  the 
factor  for  r cardation  due  to  applied  rudder  not  made  to  be  as  sever?  as 
that  indicated  by  Figure  1.  The  author  believes  that  HMAS  KIMBLA-  could  hold 
a beam  course  in  winds  above  50  kn,  but  much  depends  on  the  build-up  of  sea 
state.  HMAS  KIMBLA's  rudder  may  also  be  proportionately  larger  than  that  of 
the  ship  whose  characteristics  ere  illustrated  in  Figure  1.  Until  actual 
data  are  available,  there  is  little  point  in  pursuing  this  in  detail.  From 
the  viewpoint  of  our  subsequent  use  of  the  model  (and  Figure  8),  that  is  to 
Improve  the  dead  reckoning  in  normal  weather,  our  main  interest  lies  in  ti.e 
behaviour  in  winds  less  than  30  kn.  In  this  area  little  comparison  with 
Figure  1 can  be  made.  Our  model  Indicates,  however,  chat  HMAS  KIMBLA  should 
be  less  affected  by  light  to  moderate  winds  and  seas  astern  than  the  larger 
vessel. 

A model  that  describes  the  variations  in  HMAS  KIMBLA's  speed  threugh  the 
water  as  a function  of  weather  has  been  described.  The  only  variables  required 
as  ‘ .'ut  to  the  calculations  are  those  routinely  recorded  by  the  ship’s  officers 
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11.  RECOMMENDATIONS 

It  should  be  re-emphasised  that  this  work  was  done  in  the  manner 
described  in  order  to  make  the  best  use  of  the  available  data  from  several 
past  trials.  Improved  dead  reckoning  was  needed  in  order  to  deduce  surface 
currents.  The  obvious  way  to  improve  our  understanding  of  how  HMAS  KIMBLA 
or  any  other  ship  performs  as  a function  of  sea-state  is  to  obtain  some  direct 
measurements  of  speed  through  the  water,  along  with  a continuous  record  of  the 
pertinent  weather  data.  Thera  is  little  point  in  attempting  to  refine  this 
particular  model  without  such  data.  Thus  it  is  essential  that  HMAS  KIMBLA 

be  fitted  with  a speed  log,  as  is  the  case  for  most  ships.  If  a dual  axis 
log  is  fitted  actual  data  on  the  leeway  could  also  be  obtained.  The  weather 
data  could  be  improved  simply  by  recording  the  observations  at  the  time  of  any 
significant  change,  in  addition  to  the  regular  A hr  entry.  A continuous  record 
of  relative  wind  velocity  would  be  most  helpful.  The  main  lack  in  sea-state 
observations  concern^  data  on  wave-period  or  wave-length.  A reliable  ship- 
borne  wave  and  swell  recorder  would  be  of  value.  However,  a simpler  approach 
may  be  quite  useful  for  these  purposes.  The  major  factor  pertaining  to  loss 
of  speed  due  to  sea- state  is  probably  the  pitching  motion  of  the  ship.  A 
three-axis  accelerometer  installed  in  the  bows  of  the  ship  could  record  pitch, 
roll/yaw  and  longitudinal  decceleration  and  would  have  the  advantage  cf 
automatically  recording  the  resonance  effect.*?  associated  with  waves  and  swell 
of  certain  relative  frequencies.  These  observations  could  be  correlated  with 
measured  speed  and  a new  model  derived,  perhaps  through  the  application  of 
multivariate  analysi.s.  Such  a model  would  be  useful  to  reduce  dead  reckoning 
data  on  occasions  when  the  speed  log  wa.s  inoperative,  as  sometimes  occurs. 

Associated  analysis  would  also  be  useful  in  evaluating  and  correcting  for  errors 
in  the  recorded  log  speed  under  certain  conditions. 
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